Engineering two-mode entangled states between two superconducting resonators by 
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We present an experimental feasible scheme to synthesize two-mode continuous- variable entangled 
states of two superconducting resonators that are interconnected by two gap-tunable superconduct- 
ing qubits. We show that, with each artificial atom suitably driven by a bichromatic microwave field 
to induce sidebands in the qubit-resonator coupling, the stationary state of the photon fields in the 
two resonators can be cooled and steered into a two-mode squeezed vacuum state via a dissipative 
quantum dynamical process, while the superconducting qubits remain in their ground states. In 
this scheme the qubit decay plays a positive role and can help drive the system to the target state, 
which thus converts a detrimental source of noise into a resource. 
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I. INTRODUCTION 

Solid state devices based on superconducting cir- 
cuits are promising candidates for studying fundamental 
physics and implementing quantum information proto- 
cols These devices, which can be designed and 
fabricated on demand, provide an unprecedented level of 
tunability and flexibility in the implementation of scal- 
able quantum information processing. Three basic types 
of superconducting qubits, including charge @, ffuxQ 
and phase qubits [8(, have been theoretically and exper- 
imentally investigated to realize quantum computation 
and communication. In general, these qubits are cou- 
pled either through circuit elements such as cap acitances 
and inductances, or via a quantum bus 0, Il0j |. a super- 
conducting resonator. However, recent researches have 
shown that superconducting resonators can also be ad- 
dressed and manipulated for quantum state generation 
fllT - [l8l ] . Particularly, plenty of works have focused on the 
extension from single to more versatile multi-resonator 
architectures |19l - t23 |. where the resonators are intercon- 
nected by superconducting qubits. These multi-resonator 
systems allow manipulation of spatially separated photon 
modes and engineeri ng qua ntum states between physi- 
cally distant cavities [23-35]. 

In this work, we introduce an efficient scheme for 
the preparation of two-mode entangled states between 
two spatially separated superconducting resonators that 
are interconnected by two gap-tunable superconducting 
qubits. This proposal exploits a bichromatic microwave 
driving of each qubit to modify its coupling to the res- 
onator modes, inducing sidebands in the qubit-resonator 
coupling [3(1 H3] , and actively utilizes the qubit decay to 
drive the system to the desired state. We show that, com- 
bining the sidebands with the qubit decay, the station- 
ary state of the photon fields in the two resonators can 



be steered into a two-mode squeezed vacuum state via 
a dissipative quantum dynamical process, while the su- 
perconducting qubits remain in their ground states. We 
detail the actual physical implementation of this scheme 
with the superconducting LC resonators inductively cou- 
pled to two gap-tunable a-loop flux qubits. This work 
may open up promising perspectives for quantum infor- 
mation processing with coupled networks of supercon- 
ducting resonators. 



II. THE MODEL 
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FIG. 1. (Color online) The schematic for generating entan- 
glement between two superconducting resonators. These res- 
onators are interconnected by two gap-tunable superconduct- 
ing qubits, which are independently controlled by external 
circuits. 



We consider two gap-tunable superconducting qubits 
with ground state \g) and excited state \e), coupling to 
two superconducting resonators, as shown in Fig. 1. The 
energy levels of each qubit is individually driving by a 
bichromatic microwave field to induce sidebands in the 
qubit-cavity coupling. A typical gap-tunable supercon- 
ducting qubit is the modified flux qubit where two junc- 
tions forms a superconducting quantum interference de- 
vice (SQUID) to allow the modulation of the qubit gap 
38] . This modified flux qubit allows a fast change of the 
qubit resonance frequency while remaining at the degen- 
eracy point. The free Hamiltonian describing the two 
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qubits and the resonator modes is (let H = 1) 

2 c 
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where <5 is the qubit gap without external driving, cr A is 
the Pauli operator, za\ is the resonator frequency, and a\ 
is the creation operator for the Xth resonator. The bare 
qubit-rcsonator interaction is given by 
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where <7j is the coupling strength between the qubits 
and the resonators, and o± the qubit ladder operators. 
We now consider a a z driving with frequencies for 
the Xth qubit (37j . For flux qubits, this driving can be 
implemented by an external magnetic field that induces 
a flux driving [38| . Expressing the driving amplitude by 

the ratios f;^ to the driving frequencies, we have 
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Keeping only the leading order in the rations £ ; 
have 
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By setting u^ 1 ] 



h.c)]} + H.c. 
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V\ , u d 2 = 5 — v-i , and neglecting the fast oscillating terms 
under the rotating-wave approximation, we can get 

+H.c. (8) 

This Hamiltonian describes the red sideband (u>^\ = 

5 — v\) and blue sideband {ur d2 = ^ + ^2) excitations 
of the first qubit induced by the bichromatic microwave 



field with frequencies u d X \ and 
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and the blue side- 
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band (uJ d { = 6 + fi) and red sideband (uJ d 2 = S — V2) ex- 
citations of the second qubit induced by another bichro- 

(2) (2) 

matic microwave field with frequencies uj d x and ur d ^ . The 
schematic diagram of these engineered qubit-photon cou- 
plings is shown in Fig. 2. 



In the following, the system-environment interaction is 
assumed Markovian, and then is described by a master 
equation in Lindblad form. Here we consider only relax- 
ation and dephasing of the superconducting qubits, with 
the relaxation rate T r and dephasing rate T v . The time 
evolution of the density operator p for the whole system 
is then described by the master equation 



^ = -i[H, p] + T r D[al]p + T r D[a 2 _]p 
+r v D[al]p + T v D{a 2 z ]p, 



(4) 



where = Ho + i^ + ii/, and D[A]p = 2ApA^ -A^Ap- 
pA^A. 

We proceed to perform the unitary transformation 
Ui{t) = e- lHot , which leads to 
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To engineer the desired qubit-photon coupling, we per- 
form another unitary transformation 

U 2 (t) =Te- i ^ & dr)dr 



H d (t x )dh 



* i-ti 

dh / dt 2 H d {ti)H d {t 2 ) + ... (6) 
Jo 



FIG. 2. (Color online) The schematic of the engineered side- 
band couplings for the two superconducting qubits. The en- 
ergy gap of the jth qubit oscillates at frequencies and 
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2 induced by the flux driving of 4>^(t)(j = 1, 2). 
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g{ — 82, and 8i > 82- In this case the Hamil- 
tonian (f5J) can be rewritten as 

H = 4(8iai + 9 2 a|) + 4(8ia 2 + 82a!) + H.c. (9) 

We apply a unitary transformation p = 5^ (C)/35(C), 
where 5(C) = e (^-C&\4^ = tanlT 1 8 2 /8i. The mas- 
ter equation Q for the system now becomes 

dp 
~dt 



-i[H, p] + T r D[<rl]p + T r D[a 2 _]p 



+r !p D[*l}p + r v D[a 2 z ]p, 



where 



H=Je 2 -e 2 (a 1 a 1 + 



a2o\) 



H.c. 



(10) 
(11) 
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The master equation (flOl) is equivalent to that of two 
independent Jaynes-Cummings models of two two- level 
atoms coupled to two cavity modes respectively, with 
atomic spontaneous emission and dephasing included. In 
the transformed picture, the jth superconducting qubit 
pulls the jth resonator mode into the vacuum state by 
a process in which the qubit continuously absorbs pho- 
tons and decays back to its ground state. Therefore, the 
steady state of the system in the transformed picture 
is the vacuum |00) r for the resonator modes 01,02, and 
Iff) 1(5)2 for the superconducting qubits. Reversing the 
unitary transformation, it can readily be seen that the 
state 

m=S(Q\QQ) r \g) 1 \gh (12) 

is the unique stationary state of the master equation (U) 
i.e., 

p(t-> 00) = \ip s )(il>s\- ( 13 ) 

The steady state (fT2|) does not depend on the initial 
photon state of the resonators. Thus the resonators does 
not have to be initially cooled to the vacuum state in 
order to prepare such a state. Another distinct feature 
is that, in this scheme the qubit decay plays a positive 
role and can help drive the system to the target state, 
which thus converts a detrimental source of noise into a 
resource. This feature is particularly favorable for super- 
conducting quantum circuits, since decoherence of the su- 
perconducting qubits is one of the generic sources of noise 
and sets the limits of quantum coherence in supercon- 
ducting quantum circuits. The state e'* aia2 ~'' a i a 2 |00) r is 
a two-mode squeezed state of the photon fields in the 
two resonators, which exhibits Einstein-Podolsky-Rosen 
(EPR) entanglement. The degree of squeezing (squeeze 
parameter £) is determined by the ratio of 82 to Gi, 
which can be controlled on demand through tuning the 
external experimental parameters. 

To check the above model, we proceed to numerically 
solve the master equation including the decay of the res- 
onator modes, i.e., 

^ = -i[H, p] + T r D[aL]p + T r D[a 2 _]p + T v D[a\]p 

+T v D[a 2 ]p + mDla^p + n 2 D[a 2 ]p, (14) 

where ml\ is the photon loss rate for the Ath resonator. 
To quantify the validity of the proposal, we exploit the 
total variance V = ((Au) 2 + (Aw) 2 ) of a pair of EPR-likc 
operators u = X\ + X 2 , and v = Pi — P2 , with X\ = 
(ax+a[)/V2, and P A = -i(a\-a\)/V2,\= 1,2. A two- 
mode Gaussian state is entangled if and only if V < 2. 
For an ideal two-mode squeezed vacuum state, the total 
variance V — ((Aw) 2 + (Ati) 2 ) = 2e -2 ^, implying this 
state exhibits EPR entanglement. 

In Figure 3 we display the numerical results for time 
evolution of the total variance V together with the re- 
sults for an ideal two-mode squeezed vacuum state, un- 
der different values for the decay rates T, and k, where 



r r ~ T v ~ r, and k± ~ k 2 ~ K are assumed. 
The relevant parameters are chosen such that they are 
within the parameter range for which this scheme is 
valid and are accessible with present-day experimental 
setups. The initial state of the system is chosen as 
p(Q) — |00) r - 1 (7) 1 1 (7) 2 • From the simulation one can read- 
ily see that the model performs very well with the cho- 
sen parameters. At steady state ideal EPR entanglement 
(62/61 = 0.75, V = 0.286) between the photon fields in 
the two resonators has been established if the resonator 
decay rate satisfies k < {61, 62, T}. However, when the 
resonator decay is comparable to or larger than T, the 
ideal EPR entanglement is spoiled by resonator decay. 
Therefore, implementing this proposal with high fidelity 
requires that k < {61, 62, T}. 
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FIG. 3. (Color online) Time evolution of the total variance 
V through numerically solving the master equation (14), to- 
gether with the result for an ideal two-mode squeezed vacuum 
state. The relevant parameters are chosen as £^ ~ ft( 2 ' — 
ft = 0.2, f« * tf> - ft = 0.15, g^ * ff W * * = 
9- 



III. TWO-MODE ENTANGLEMENT BETWEEN 
TWO SUPERCONDUCTING LC RESONATORS 

In principle, the generic model can be implemented 
with various superconducting resonators coupled by two 
Josephson-junction-based superconducting qubits. Here 
we propose a convenient demonstration with two su- 
perconducting LC resonators inductively coupled to two 
gap-tunable flux qubits, as shown in Fig. 4. 

The superconducting LC resonators made of a capaci- 
tor and an inductor are described by a simple harmonic 
oscillator Hamiltonian [3^, H(| H r = 5Zx=i ^(0^6^ + 
1/2), where the resonance frequency v\ is determined by 
the respective capacitance C\ and the inductance L\: 
V\ = \j\JL\C\. For a micrometer scale LC resonator, 
the resonance frequency v\ is on the order of several 
GHz[39ll40T|. The quality factor Q for a LC resonator can 
reach Q > 10 3 , which leads to a decay rate for the res- 
onator on the order of several MHz. These LC resonators 
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are interconnected by two gap-tunable flux qubits which 
are inductively coupled to the resonators. 



of the qubits, {\g), |e)}, the total Hamiltonian reads 
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FIG. 4. (Color online) The schematic of two superconducting 
LC resonators inductively coupled to the a-loop flux qubits. 
The Ath flux qubit, controlled by the static magnetic flux & x ' 
and time-dependent flux driving (j>^(t), is coupled to the ith 
resonator via mutual inductance M^ x \ 

We consider the so-called tunable-a flux qubits, as 
demonstrated in Ref. [H|. Different from the flux qubit 
with three Josephson junctions, one of which has the cou- 
pling energy smaller than that of the other two junctions 
by a factor a, here the small a junction is replaced by 
a so-called a loop, formed by a SQUID with two iden- 
tical Josephson junctions. The Ath flux qubit can be 
operated at the degeneracy point with an external ap- 
plied magnetic flux of $' A ^ ~ 3>o/2 and behaves effec- 
tively as a two- level system, where $0 = h/2e is the flux 
quantum. The qubit gap can be controlled via the flux 
driving (t) through the SQUID loop, which can be 
implemented by the external microwave control lines. In 
this case, the Hamiltonian for the Ath flux qubit can be 
written as 



#(A) 



1 



{eaM+A[<t>M(t)]*W} (15) 



in the basis of the persistent current states {| O), | O)}. 
Here e = 2/ p (<I>( A ) — $o/2), with I p is the persistent cur- 
rent in the qubit loop. The qubit gap A[<// A ) (t)} de- 
pends on the flux driving, and can be separated into 
a static part S and a time-dependent oscillating part 
cx cos(uj^t). 

The two superconducting LC resonators couple to the 
flux qubits via the mutual inductance. The interaction 
between the Ath qubit and the Ith LC resonator can be 
described by a coupling of dipolar nature in the basis of 
the persistent current states 



rrOO _ .(^) (A) / - 1 , A \ 
H I ~ 9l a z [ a t + a 0: 



(16) 



where the strength of the coupling is g^ — M^Ip^J-^, 

with the mutual inductance between the Ath qubit 

and the ith LC resonator. In the basis of the eigenstates 



H = 2_y77<J z +v x a , x a\ 



At; 



A=l 



]T 5i (A) (cos0a A -sm6a x x ){a\ + o,)] (17) 



1=1 



with uj q = Ve 2 + A 2 and tan 9 = A/e. At the degeneracy 
point e ~ 0, this leads to 



A=l 



1=1 



^ 5i (A) (a A +a A )(a[+aO} 



(18) 



1=1 



From this Hamiltonian, one can engineer the desired 
photon-qubit couplings and prepare the photon fields in 
the two LC resonators into the two-mode squeezed vac- 
uum state, via a dissipative quantum dynamical process 
following the reasoning in the above section. 

We now discuss the relevant experimental parameters. 
Taking the capacitance and inductance of the LC res- 
onators as 12 pF and 250 pH leads to vi ~ 2/2 — 2.9 
GHz [4lJ]. The mutual inductance between the Ath 

flux qubit and the Ith LC resonator can be about 20 pH 
[4lj ] . Thus we can obtain the coupling strength g ~ 200 
MHz. If we take £1 = 0.2,^2 = 0.15, then we obtain the 
effective coupling strength 81 ~ 40 MHz, and 02 — 30 
MHz. The time for preparing the steady state (|12p is de- 
termined by the decay rate T of the flux qubit. Provided 
that yOf — ©f — r, this time will be the order of a few 
times 1/r. Moreover, if k < T, the generated state will 
be nearly an ideal two-mode squeezed vacuum state. If 
we take T ~ 20 MHz for the flux qubit, the preparing 
time will be about 250 ns. 



IV. CONCLUSION 

To conclude, we have presented an efficient scheme for 
the preparation of entangled states between two super- 
conducting resonators that are interconnected by two 
gap-tunable superconducting qubits. We have shown 
that, with each qubit individually driven by a bichro- 
matic microwave field to induce sidebands in the qubit- 
cavity coupling, the stationary state of the photon fields 
in the two resonators can be steered into a two-mode 
squeezed vacuum state via a dissipative quantum dynam- 
ical process. This proposal actively exploits the qubit 
decay to drive the system to the desired state and does 
not depend on the initial photon state of the resonators, 
which can be implemented with superconducting LC res- 
onators inductively coupled to a-loop flux qubits. 
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